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We calculate the ground-state two-spin correlation functions of spin-i quantum Heisenberg chains 
with random exchange couplings using the real-space renormalization group scheme. We extend the 
conventional scheme to take account of the contribution of local higher multiplet excitations in each 
decimation step. This extended scheme can provide highly accurate numerical data for large systems. 
The random average of staggered spin correlations of the chains with random antiferromagnetic (AF) 
couplings shows algebraic decay like 1/r 2 , which verifies the Fisher's analytic results. For chains 
with random ferromagnetic (FM) and AF couplings, the random average of generalized staggered 
correlations is found to decay more slowly than a power-law, in the form close to l/ln(r). The 
difference between the distribution functions of the spin correlations of the random AF chains and 
of the random FM-AF chains is also discussed. 



I. INTRODUCTION 

One-dimensional (ID) quantum spin systems have attracted much attention over the decades. This is not only 
because these systems have been a good testing ground for various theoretical techniques and approximations but 
also because they exhibit a wealth of fascinating phenomena in their ground states and low-lying excitations. These 
include quasi long-range-order (LRO), topological order and ground-state phase transitions, which are all purely 
quantum effects due to the low-dimensionality of the systems. Among these quantum phenomena, the effects of 
randomness on quantum spin systems have been studied intensively by many groups. These studies revealed the 
appearance of various exotic phases, which are realized neither in regular quantum systems nor classical random 
systems. The interplay of randomness and quantum fluctuations plays an essential role jh. these phases. Systems 
which are gapless in the absence of randomness are unstable against weak randomness,IM [While gapped systems, 
e.g. integer spin chains in the Haldane phase and dimerized chains, are comparatively robust. uu 

A simplest model of the ID random spin-i Heisenberg spin systems is described by the Hamiltonian of the form 



L-l 

h = Jis t ■ s* i+ i, (i) 



where s$ are S = h spin operators and the exchange coupling constants Ji are distributed randomly according to 
a probability distribution P(Ji). There are several quasi-lD systems which provide the realizations of the model 
Hamiltonian (|l|). To our knowledge, the first example of such systems belongs to the class of the organic charge- 
transfer salts tetracyanoquinodimethanide (TCNQ) .til The low-temperature magnetic properties of these systems are 
successfully described by the model Hamiltonian (Q) with random antiferromagnetic (AF) c»upling,Lj in which the 
couplings Ji are restricted to take positive random values. A more recently studied systems is SrsCuPti-^Ir^Oe- 
While the pure compounds Sr3CuPtC>6 and SraCuIrOe represent, respectively, AF and ferromagnetic (FM) spin 
chains, Sr3CuPti_ x Ir K 06 contains both AF and FM couplings. The fraction of FM bonds is simply related to the 
concentration x of Ir. These compounds are modeled adequately by the Hamiltonian (Q) with P{Ji) — pS(Ji + J) + 
(1 —p)6(Ji — J), where p is the probability of FM bonds. The experimental data of the susceptibility of the compounds 
were in fact explained-successfully by a theory based on the Hamiltonian (Q) with the above probability distribution 
of FM and AF bpnds.t] The model (|l|) is also realized in the low-temperature regime of randomly depleted ID spjn-i 
even-leg ladders.EEl In these systems, effective S = | spins are induced in the vicinity of each depleted siteJiirO If 
the density of depleted spins is sufficiently small, the induced effective spins are the only degrees of freedom which 
are active in the low-energy limit. The strength of the residual interaction between effective spins depends on the 
distance between the effective spins, and the coupling can be either antiferromagnetic or ferromagnetic, depending 
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on whether or not two spins are depleted from the same sublattice. Thus, the low-energy physics of the randomly 
depleted spin-i even-leg ladders is described as a random spin chain with both AF and FM exchange couplings. 13 

For ID random spin chains containing only AF couplings, rather complete theoretical understanding on their low- 
energy properties has been achieved. One of the most powerful techniques to study such systems is the real-space 
renormalization group (RSRG) method introduced by Ma, Dasgupta, and Hu.El The basic idea of this method is 
iterative decimation of spin degrees of freedom by integrating out the strongest bonds in the chain successively. 
For the random AF chains this procedure keeps the form of the Hamiltonian (|l|) and renormalizes the probability 
distribution P(J,). Fisher has given a solution to this RG equation of P(Ji) that becomes exact in the low-energ, 
limit and shown that any normalizable initial distribution P{Ji) flows to a single universal fixed point distribution 
The ground state of the chain which belongs to the fixed point is characterized by the "random singlet phase," where 
each spin forms a singlet with another spin which can be arbitrarily far away. From the intuitive picture of the random 
singlet phase, Fisher has also shown that the random average of the two-spin correlation function in this phase decays 
algebraically like 1/r 2 with r being the distance between the two spins. The main contribution to the average comes 
from the rare events that two spins separated by the distance r form a singlet pair. The 1/r 2 power-law decay has 
been verified by a numerical calculation for the random XX spin-Jj chains,Ea which can be mapped to a disordered 
system of noninteracting fermions. 

The random spin chains containing both AF and FM couplings have also been studied theoretically for the last 
five years, from which a qualitative picture of the low-energy physics of such random FM-AF chains has emerged. 
Westerberg et al£3 have adapted the RSRG scheme and shown through extensive numerical simulations that the 
distribution P{Ji) is renormalizcd to a single universal fixed-point distribution unless the initial distribution is highly 
singular around Ji = 0. The chain at this fixed point can be viewed as an ensemble of weakly interacting large 
effective spins whose average size S oc T~ K with n » 0.22 for temperature T — > 0. These large effective spins 
are generated as a result of decimations of two spins coupled via a strong FM coupling into, a larger spin. The 
results of the RSRG simulations are supported by a recent calculation by Frischmuth et al.^B who have used the 
continuous-time quantum Monte Carlo loop algorithm. In contrast to the success in the quantitative calculations on 
the thermodynamic properties, less is known about the spin correlations in the random FM-AF chains. Since the 
average size of the effective spins becomes very large in the low-energy (long-distance) limit, one can expect that the 
system should be close to a classical spin chain that can show the LRO of the generalized staggered spin component. 
One might also argue, however, that there should be no LRO in ID quantum spin chains. It is therefore an interesting 
open question how the correlation function of the generalized staggered moment (whose definition will be given in 
Sec. Ill) behaves in the ground state. The purpose of this paper is in fact to present results of our extensive numerical 
calculation of the two-spin correlation function. We find that it decays very slowly with the form close to l/ln(r), 
which is consistent with the naive argument that the ground state is extremely long-range correlated, but not really 
long-range ordered. 

The outline of this paper is as follows. Section II A is devoted to a brief review of the RSRG scheme of the generalized 
version which is applicable to both random AF and FM-AF case. Using the Wigner-Eckart theorem,! 2 -!! we si mplify 
the algorithm to allow for calculating the two-spin correlation function between original S — \ spins. In section II B , 
in order to achieve higher accuracy, we extended the RSRG scheme to take into account the contributions of local 
excitations to the ground state of the whole system. We perform these conventional and extended RSRG algorithm 
numerically on both random AF and FM-AF chains to calculate the "staggered" spin correlations on their ground 
states. The results on the random AF case are shown in section III A . In section fill B , we show the results on the 
random FM-AF case. Analyzing the data obtained with the extended RSRG method (which can be applied to larger 
systems than the DMRG), we conclude that the mean correlations on the random FM-AF case decay very slowly with 
logarithmi c dep endence on r. We also discuss the distribution functions of the logarithm of the correlation functions 
in section [II C. We find that, in the random AF case, the rare spin pairs, which are strongly correlated, dominate 
the mean correlation while such rare events are not essential in the random FM-AF case. Finally, our results are 
summarized in section IV. 



II. THE RSRG ALGORITHM 



Ma, Dasgupta and Hu introduced a RSRG procedure to investigate the low-temperature properties of random AF 
spin chains.El The method has been generalized to the random FM-AF case by Westerberg et alJtB In this section we 
explain our extension of this scheme to calculate the ground-state two-spin correlation functions. We begin with a 
brief review of the RSRG method for the random FM-AF case. 
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A. Conventional RSRG 



Let us consider a random FM-AF spin-i chain described by the Hamiltonian (|l|). The basic strategy of RSRG is to 
decimate spin degrees of freedom by combining two spins connected via a strong bond into one effective spin. Conse- 
quently, the system is described in terms of effective spins of various sizes coupled by random exchange interactions, 
although the original Hamiltonian ([j]) consists of only S = 1/2 spins. We accordingly treat the effective Hamiltonian, 

H = Y / JiSi-S l+1 , (2) 
i 

where both the coupling Ji and the size of the effective spins Si are random. We call the S =1/2 spins Si appearing 
in the Hamiltonian ([!]) "original" spins in the following to distinguish them from the effective spins Si. 

We define A; as the energy gap between the ground-state multiplet and the first excited multiplet of the corre- 
sponding bond Hamiltonian Hi — JiSi ■ Si+i, 

_ r \Ji\(s t + s l+1 ) (Ji<o) 

^ \Ji(\Si- S i+1 \ + l) (Ji>0) 

and focus on the bond with the largest gap A; in the chain. The terms in the effective Hamiltonian (||) which involve 
the effective spins Si and Si+i arc 



where 



H' = H' a + H[ (3) 

H' = JiSi-Si +1 , (4) 
H[ = Ji-iSi-i ■ Si + Ji + \Si + i ■ Si + 2- (5) 

If Ai is much larger than the gaps of the neighboring bonds, A;_i and A; + i, the spins Si and to a good 

approximation, are frozen into the ground-state multiplet of the local Hamiltonian H . We. therefore, replace the 
block composed of Si and Si+i by the single effective spin S. The Wigner-Eckart theorenO then implies that both 
Si and Sj+i are proportional to S: 

Si = aS, 

Si +1 = PS, (6) 

where a and (3 can be obtained from the Clebsch-Gordan coefficients. Substituting Eq. (^) into the four-spin Hamil- 
tonian (Q), we obtain, apart from a constant term coming from Hq, 

H = Ji-iSi-i ■ S + Ji+iS ■ Si + 2, (7) 

where 

J;-i = aJi-i, 
Ji+i = P'h+i- 

The case where Ji is antiferromagnetic with Si = Si+i needs a special treatment. In this case the two spins Si and 
Si+i form a singlet, and accordingly, we remove both spins from the effective Hamiltonian. Between the spins Si-i 
and Si + 2, a coupling is generated through a second-order process that virtually breaks the singlet of Si and Si + i. We 
obtain 

H = JSi-i ■ S1+2, (8) 

where 

oJi 
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By replacing the four-spin Hamiltonian (H) with H [Eq. (f7j) or Eq. (||)] in the Hamiltonian of the whole system, 
we obtain a new effective Hamiltonian (See Fig. [|). We note that this procedure preserves the form of the effective 
Hamiltonian (0) but changes the distributions of the exchange couplings, J/, and the size of effective spins, Si. We 
repeat this procedure of integrating out the strongest bonds in a chain successively until the distribution functions for 
J\ and Si converge to a scaling form._T-his RG flow has been investigated extensively for various initial distributions 
including both AF and FM-AF case.Elll3 In particular, for the random AF S = 1/2 chains Fisher has solved the RG 
equation exactly.El 

It is also possible to calculate the correlation functions between original spins within the scheme described above. 
Here we make use of the fact that the original spin operator is, at every step in the RSRG procedure, proportional 
to the effective spin Si to which it belongs. We can keep track of the coefficient for each original spin operator by 
applying Eq. (Q) at each step. At the step where the effective spins Si and Si+i is added into 5, we calculate the 
correlations between si and Sj in the ground state of the bond Hamiltonian H' a , 

(s*j • 8j) = aiOtj(Si ■ Si+i), (9) 

where (• • ■) represents the expectation value in the ground state; Sj and Sj belong to Si and Sz+i, respectively; OLi and 
Oj are the proportionality coefficients for each spin. 



B. Extension of the RSRG scheme 



As seen in the previous subsection, the conventional RSRG procedure consists of "diagonalizing the bond Hamil- 
tonian with the largest gap" and "projecting the low-energy states onto the lowest multiplct." This means that we 
completely neglect the contribution of the excited multiplets of the local Hamiltonian (3) to the ground state of the 
whole system. This approximation is valid only if the energy gap of H' is much larger than the ones of the neighbor- 
ing bonds. Fisher's solution for the random AF S = 1/2 chain becomes asymptotically exact since this condition is 
satisfied near the fixed point of the RG flow. 

However, the condition is often not satisfied, in particular, in the early stage of the RG, unless the initial distribution 
of energy gap is very broad. As a result, we have no choice but to cut off the contributions of local excitations. This 
is a poor approximation that has a serious effect especially on the calculation of the expectation values of microscopic 
operators such as the original spins. The ground-state correlation functions between original spins calculated via the 
conventional RSRG in fact deviate largely from those obtained from the DMRG method and the extended RSRG 
algorithm described below (See Fig. [| and Fig. ||in section III). 



One possible prescription to avoid this error is "to keep the local multiplet excitations" at each step. Of course, if 
we keep all eigenstates of the bond Hamiltonian Hq at every step of the RG, the calculation on H' in the final step 
is equivalent to the exact diagonalization on the Hamiltonian of the whole system. In practice what we should do is 
to extend the RSRG algorithm under the policy that we keep as many states as we can store in computer memories. 
In the original RSRG scheme a segment of original spins are combined and represented by a single effective spin. 
In the extended scheme we keep more states than the lowest multiplet and call the segment a 'block.' Each block 
consists of several original spins and is represented by 'block states,' which are the m-lowest eigenstates of the block 
Hamiltonian. Since we have to keep or discard all states of a multiplet to ensure the SU(2) symmetry, the actual 
number of kept states for block I is rn ; * < m. An original spin operator in block I is represented by m ; * x m ; * matrix 
on the set of the block states accordingly. 

Let us consider the effective Hamiltonian 

H = Y J H\ B) +Y J Hi,i +1 (fO) 
i i 

h ij+1 = Jis[ R) -s[+l (n) 

where H\ is a block Hamiltonian of the Zth block, diagonal in the block states; -H^j+i is a coupling Hamiltonian 

between the Zth and (Z+l)th blocks; and sf are original spin operators on the right and left edge of the Zth block, 

respectively; Ji is a coupling between Sj and sfS- In the extended RSRG scheme we renormalize the two- block 
Hamiltonian 

H Li+i = H \ B) + H i,i+i + H i+l I 12 ) 

with the largest gap A; into one block Hamiltonian (see Fig. |J). Here we define A; as the energy difference between 
the highest energy in the eigen-multiplets of H[ [+1 which will be kept and the lowest energy in the multiplets which 
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will be discarded after the decimation. The basic scheme of the extended RSRG is the same as the conventional 
RSRG except the changes described above. The algorithm is summarized as follows: 

(i) Focus on the bond with the largest gap A;. Construct the two-block Hamiltonian (|l2|) of the bond. 

(ii) Diagonalize the two-block Hamiltonian to find a set of eigenvalues and eigenstates. At this point, we can 
calculate expectation values of various operators in the two blocks, such as the two-spin correlation functions 
between the original spins, using the ground state of the two-block Hamiltonian. 

(iii) Discard all but the lowest m*(< m) eigenstates in the block Hamiltonian. 

(iv) Express operators, such as the block Hamiltonian itself and the original spin operators in the new block, in 
terms of the new block states. 

(v) Rewrite the coupling Hamiltonians between the new block and its neighboring blocks in terms of the new s 1 -^ 
and s^ R \ Diagonalize them to update the distribution of the energy gap A. 

(vi) Return to the step (i). 

We obtain all two-spin correlation functions by repeating this procedure until the whole chain is finally represented 
by one block. 



III. NUMERICAL RESULTS 



Using both the conventional RSRG and the extended RSRG with various values of m, we calculated the correlation 
functions (s*j • Sj) in the ground state of open chains for both the random AF case and the random FM-AF case. 
The maximum size of the chains used for the conventional and the extended RSRG simulations is L = 100000 and 
L = 1000, respectively. In the FM-AF case, the random average of the spin correlation, (s*j • Sj), always decays 
exponentially^ and is not an appropriate quantity to characterize the spin correlation because the sign of (s, • Sj) can 
be either positive or negative depending on the number of AF bonds between Si and Sj. Instead we introduce the 
"generalized staggered" correlation function 

CQi-j\) = Vij ft-3 j ), (13) 

where r/ij = Ylk=i s § n (~Jk) f° r 3 > For random AF case, C(r) is the usual staggered spin correlation. We 
take the random average of C(r) and lnC(r), which represents the mean correlations and the logarithm of the 
typical correlations, respectively. Note that it is impossible to take the random average of lnC(r) numerically in the 
conventional RSRG algorithm, within which C(r) =0 for two spins that do not form a singlet pair. To check the 
results pithe RSRG methods, we also calculated the correlation functions onL= 100 open chains using the DMRG 
methodE3 with the improved algorithm proposed by White.Ej The number of kept states in the DMRG calculation was 
up to 100 and 200 for the random AF and FM-AF case, respectively. In both cases, the mean and typical correlations 
calculated by the extended RSRG are in excellent agreement with those by DMRG (see below). We note that the 
systems we have treated are much bigger than those in the earlier work by Hida,u in which the DMRG method is 
applied for the FM-AF chains. 



A. Random AF case 



As a typical probability distribution of the random AF case, we choose the box- type bond distribution P(Ji), 

P ( J *) = {o° (otherwise) ^ (14) 

where the cutoff, Jq, is taken as energy unit. For the random AF case, it is known that any normalizable initial 
distribution flows to a single universal fixed point.El Hence, the results we obtained for the initial distribution ( |l4| ) are 
generic. The number of sample chains we used for each method are shown in Table [|. 

Before discussing our numerical results, we briefly comment on the m-dependence of the data of the extended 
RSRG. As noted in the last section, the extended RSRG can provide more accurate results as the number of kept 
states, m, increases. In the random AF case, the ground-state multiplet of a block Hamiltonian is either singlet or 
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doublet, depending on the number of original spins belonging to the block. The degeneracy of each low- lying excited 
multiplet is, therefore, always small. As a result, we can keep considerably many multiplets even if m is rather small. 
We estimate the m-dependence of the results of the extended RSRG from the numerical data with m — 10, 20, 30. 

Figures || and || show, respectively, the numerical data of the mean correlations, ((C(r))), and the average of the 
logarithm of correlations, ((lnC(r))), where ((• • •)) represents random average. It is clear that the data of the extended 
RSRG converge already at m = 30, and we can regard those data with m = 30 as those in the limit m — > oo. In 
Fig. H the data of the extended RSRG with m = 30 are in good agreement with those of DMRG for r < 30, where 
the DMRG data are considered to be exact and free from finite-size effects. On the other hand, the conventional 
RSRG largely underestimates the correlations, but its data are on a line parallel to the data of the extended RSRG 
in the log- log plot. We conclude from these observations that the results of the extended RSRG are quantitatively 
reliable, whereas the conventional RSRG can be used to estimate the exponent of the power-law decay. Encouraged 
by the agreement between the extended RSRG and the DMRG data, we anticipate that the extended RSRG provide 
quantitatively reliable data even for r > 30 where the reliable DMRG data are not available. From the data of the 
extended RSRG with m = 30 for r < 300, where the data are expected not to be hampered by finite-size effects, we 
estimate the asymptotic form of the average correlations to be 

«C(r)))~r- 2 . (15) 

For the average of In C(r), we also rely on the data obtained from the extended RSRG scheme. Figure || gives 

<(lnC(r)))~-r - 5 . (16) 

Both results (|l5| ) and (|l6|) agree with Fisher's theoryU and can be considered as a numerical-aerification of his solution 
on the mean and typical correlations for th£_Heisenberg case. For the random XX chainail and for a related model 
of the random transverse-field Ising model,E£l the power-law behavior (|l5|) and (^) is observed numerically. To our 
knowledge, Figs. and are first numerical results confirming Fisher's theory for the Heisenberg case. 



B. Random FM-AF case 

Westerberg et a/.El showed that the RG trajectories of the random FM-AF chains flow towards a single universal 
fixed point under the conventional RSRG procedure, unless the initial distribution of couplings is more singular than 
P(Ji) ~ |Ji| _J,c , y c ~ 0.7. In this section we investigate the spin correlations at this fixed point with the extended 
RSRG scheme. For this purpose we assume the box-type bond distribution P(Ji), 

[0 (otherwise) , 

with Jo as energy unit, as a representative of distributions with no singularity at Jj = 0. We expect that our results 
obtained for the initial distribution ( |l7j ) should capture the universal behavior for the random FM-AF chains which 



are in the basin of the universal fixed point found in Ref. 19. We have numerically calculated spin correlations using 



the conventional RSRG, the extended RSRG with m = 30, 40, 50, 60 and the DMRG method. 

In the random FM-AF chains the degeneracy of the lowest multiplet of a block becomes larger on average as the size 
of the block grows. At the point when the degeneracy exceeds the number of kept states, m, determined beforehand 
in the algorithm, the extended RSRG scheme breaks down because we need to keep all the degenerate states in the 
lowest multiplet to preserve the SU(2) spin symmetry. In fact, we could complete the extended RSRG procedure 
without exceeding the limit of the number of kept states m — 30 (m = 60) for 30% (75%) of the sample chains 
(L = 1000). We then used only those samples for which the RSRG could be completed to take the random average. 
Although this sorting out of sample chains may lead to a systematic underestimate on the average values, we believe 
that we can correct it by carefully checking the m-dependence of the data. The number of sample chains we used to 
take the random average for each RG scheme is shown in Table |l[ 

The numerical results for the mean correlations, ((C(r)}}, are shown in Fig. || in a log-log plot. It is clear that 
the data of the extended RSRG have almost converged with m = 60 for r < 500, where the data are expected to be 
free from the effect of the open boundaries. The data with m = 60 are also in excellent agreement with the data of 
DMRG. We, therefore, regard the results of the mean correlation with m = 60 as essentially converged. We notice 
that the curve in the log- log plot are bent upward, indicating that the mean correlations ((C(r))) decay more slowly 
than a power-law. To analyze the r-dependence of the mean correlations, we plot the inverse of ((C(r)}} as a function 
of lnr in the inset of Fig. We find that the data lie on a straight line in this plot, from which we speculate that 
the mean correlations decay with the logarithmic form, 
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((c(r))) ~ r^r-y 

ax{r/r ) 



(18) 



where-fl and r$ are constants. We note, however, that Eq. ( jig ) is not the only form that can account for the numerical 
dataEj In any case the mean correlations show very weak r-dependence, probably through the form of ln(r/ro), 
certainly different from power-law. 

Figure H shows the numerical results of ((In C{r))). The data of the extended RSRG with m = 60 exhibit a similar 
behavior to the logarithm of the mean correlations; The typical correlation exp [((lnC(r)))] decays again more slowly 
than a power-law. This leads us to plot the inverse of the typical correlations as a function of lnr (see the inset of 
Fig. |^). Although the curve at m = 60 seems almost linear for 50 < r < 400, we cannot determine the r-dependence 
of the typical correlations from the figure due to the slow convergence of the data with increasing m. This slow 
convergence arises from the fact that the numerical estimate of In C(r) is very sensitive to small fluctuations of C(r) 
especially when the value of C (r) is extremely small. Calculations with much larger values of m would be needed to 
obtain the data accurate enough to determine the r-dependence of ((In C(r))). 



C. Distribution of the correlation functions 



In order to make a distinction between characteristics of ground-state correlation functions in the random AF chains 
and in the random FM-AF chains- we analyze the distributions of the logarithm of the correlations, D(x;r), where 
x = In C(r). Henelius and GirviriiS have shown numerically that for large r the distribution function of XX chains 
with random AF couplings scales to a fixed-point distribution of the form 

D(x;r) = f(r)F(x/g(r)) (19) 

with 

f(r)g(r) = 1, (20) 

g(r)=c\((lnC(r)))\, (21) 

where c is a positive constant. The scaling function F(x/g(r)) satisfies the normalization condition J F(y)dy = 1. We 
will demonstrate that the distribution D(x; r) in random Heisenberg chains also exhibits the scaling behavior Eq. ( |l9| ) 
for both the random AF and the random FM-AF case. 

Figure [?] shows the data of D(x; r) for the random AF case obtained using the extended RSRG with m — 30. 
According to Eqs. ( [l6| ) and (pl|), we can set g(r) — r 5 and f(r) — r °rn The data points (circles and squares) in 
Fig. []] collapse on a single curve, indicating that the scaling ( ^9|) applies.EZI In Fig. we also plot 

w(x,r) =wm- (22) 

which measures the contribution to the mean value of the correlation ((C(r))). Although the curves are rather 
rough due to a statistical error, it is clear that VF(x; r) has a considerable weight in the range where D(x; r) is very 
small. This means that a very few spin pairs that have much stronger correlations than typical ones give dominant 
contribution to the mean correlation ((C(r))). We also find that, as r increases, the region of x/g(r) where W(x;r) 
is large moves towards x/g(r) = 0. Indeed this is the behavior expected from Eq. (|T5|); the value of x/g(r) where 
W(x;r) takes a large weight should change as ln((C(r))) /g(r) ~ lnr/r 5 — > as r — > oo. Thus we regard the results 
shown in Fig. ^ as a further support for the random-singlet picture of the ground state of the random AF Heisenberg 
chains, where each spin forms a singlet pair with another spin that can be arbitrarily far away. The mean correlation 
((C(r))) is dominated by the contribution from the rare case in which two spins separated by distance r form a spin 
singlet. 

Figure || shows the numerical data of D(x; r) of the extended RSRG with m — 60 for the random FM-AF chains. 
Here we set g(r) — ((lnC(r)))/((ln C(r = 200))) and f(r) = l/g(r). The data points (circles and squares) for 
100 < r < 500, where we may ignore the boundary effect, clearly lie on a single scaling curve. In this range of r, 
however, g(r) changes by several percent only, from 0.95 to 1.04. To establish the scaling behavior for wide range of 
g(r), the calculations on much (exponentially) larger systems might be necessary. Such large-scale calculations are 
obviously unfeasible with computers available at present, and thus we can only conclude that our results shown in 
Fig. H are consistent with the scaling hypothesis (|l9|). In the following discussion, we regard D(x; r)//(r) in Fig. || as 
a fixed-point form of the distribution function. 
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Figure || clearly shows that the distribution function of the random FM-AF chains has a quite different form 
from that of the random AF chains, ft is essentially zero for x/g[r) > —1 and increases approximately linearly at 
x/g(r) < —1. The weight function W(x;r) representing the contribution to the mean correlation is also shown in 
Fig. |[ In contrast to the random AF case, W(x; r) has most of its weight in the region where D(x; r) is not negligible. 
This feature highlights the different nature of the spin correlations in the random FM-AF chains. As shown in the 
RSRG analysis, many spins in the random FM-AF chains correlate and form a large effective spin. This suggests that 
a spin is correlated with many other spins that belong to the same large effective spin, and, therefore, the mean value 
of the spin correlation function is not at all dominated by the "rare" events that two spins, far apart from each other, 
form a singlet pair. 

From the observation that D(x;r) j /(r) of the random FM-AF chains is negligible for x/g(r) > —A (A rs 1 in 
Fig. ^) and has an approximately linear dependence for x/g(r) < —A until it takes a maximum value, we can also 
draw a conclusion that ln((C(r)}} and ((lnC(r))) should have a similar dependence on r, in agreement with Figs. |^ 
and |^. Let us assume that the scaling function has the form 



-k(y + A) 







y< 
v > 



-A, 
-A, 



(23) 



where y = x/g(r), k and A are positive constants. Since g(r) — > oo in the limit r — > oo, the mean correlation ((C(r))) 
is dominated by the contribution from the region of small \x/g(r)\. Thus, we are allowed to use Eq. ( |23| ) for calculating 
((C(r))), even though this form is not correct for large \x/g(r)\. Using Eqs. (|l9|), ( pp| ) and (]23|), we calculate the 
mean correlation: 



«C(r)» = / e x D(x;r)dx = -k 



-Ag{r) 



g( r ) 



A] dx = 



WW 



,-Ag{r) 



yielding 



ln((C(r)» = -Ag(r) + 0( In. g(r)) = Ac((\nC(r))) + 0( In | ((In C(r)» |) (24) 

in the limit r — > oo. This result is in agreement with our observation that the curves in Figs. || and || look very similar. 
It is also consistent with the fact that W{x;r) in Fig. || stays almost in the same region of x/gir) with increasing 
r. Equation (24) is not changed qualitatively even when F(y) has an algebraic form F(y) oc [—(y + A)] a , as far as 
A > 0. Note that Eq. (E4|) does not hold in the random AF chains, for which we should take A = and a = 3. 



IV. CONCLUSION 

We have studied the spin correlations in spin-^ random Heisenberg chains for both the random AF and the random 
FM-AF case. One of the most powerful methods for the study of random spin chains is the real-space RG in which 
the spin degrees of freedom are decimated by integrating out the strongest bonds successively. In order to calculate 
the two-spin correlation functions between original S = \ spins, we modified the algorithm of the conventional RSRG 
using the Wigner-Eckart theorem. We also extended the RSRG scheme by keeping the local excited multiplets as 
block states to keep track of quantum effects. We demonstrated that the extended RSRG algorithm is very powerful, 
in the sense that it can be applied to rather large systems, and provides accurate numerical data that are in excellent 
agreement with those of the DMRG method. The numerical data of the mean and typical correlations for the random 
AF chains verified Fisher's prediction [Eqs. ([l5]) and (|l6|)] for the Heisenberg case, which cannot be mapped to 
noninteracting fermions. For the random FM-AF chains, we found that the mean correlation, ((C(r))), decays very 
slowly with the logarithmic r-dependence. We, therefore, conclude that the generalized staggered spin correlation in 
the random FM-AF chains has no LRO, although it shows very long-range correlations decaying slower than power- 
law behavior. Investigating the distribution of the logarithm of the spin correlation functions, we also clarified the 
different nature of the ground states of the random AF and FM-AF chains. In both cases, our data of distributions 
satisfy the scaling hypothesis, Eq. (jl^). Analyzing the form of the scaling function, we found that the "rare spin 
pairs", which are correlated much stronger than typical ones, dominate the mean correlation ((C(r))) in the random 
AF case. It is also shown that the r-dependence of the mean correlations in this case is quite different from that of 
the typical ones. These features strongly support the fact that the "random singlet phase" is realized in the ground 
state of the random AF Heisenberg chains. On the other hand, the scaling form of the distribution functions for the 
random FM-AF chains suggests that such "rare spin pairs" do not play an important role in these chains. We also 
deduced that ln((C(r))) and ((In C(r))) have very similar r-dependence, using the scaling hypothesis and the specific 
feature of the scaling function we obtained numerically. This result is non-trivial and consistent with our numerical 
data (Figs. | and ||). 
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RG scheme 



number of samples 



conventional RSRG 200 

extended RSRG (m = 10) 2460 

extended RSRG (m = 20) 1878 

extended RSRG (m = 30) 1899 

DMRG 1887 



TABLE I. The number of samples simulated to take the 
random averages in the random AF chains. 



RG scheme number of samples 

conventional RSRG 150 

extended RSRG (m = 30) 1319 

extended RSRG (m = 40) 1446 

extended RSRG (to = 50) 1890 

extended RSRG (to = 60) 1557 

DMRG 204 



TABLE II. The number of samples simulated to take the 
random averages in the random FM-AF chains. 
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J/-i , h J/+i 

/ v / 



J/-1 
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FIG. 1. Schematic picture of renormalization procedure in the conventional RSRG. 





FIG. 2. Schematic picture of renormalization procedure in the extended RSRG. The dashed and solid lines represent eigen- 
states of each block Hamiltonian which are discarded and kept, respectively. 
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The statistical errors of the data of the extended RSRG for r < 300 and of the conventional RSRG is smaller than the size of 
symbols. 
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FIG. 4. ({In C(r))) versus r ' 5 for the random AF case. The statistical errors of the data are smaller than the size of symbols. 
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+ : extended(m=50) 
x : extended(m=40) 
: extended(m=30) 



—i i i i i L 



—i i i i i i i l_ 



j i i i i i i i 



10 



10' 



10' 



FIG. 6. ({In C(r))) of the random FM-AF case as a function of r. The statistical errors of the data for r < 500 are smaller 
than the size of symbols. Inset: exp[— {(In C(r))}] versus r. 
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FIG. 7. Probability distributions of the correlation functions in the random AF case for r = 60,100,200. W(x;r) is also 
plotted for r = 60, 100, 200 with a solid, dotted, dashed curve, respectively. 




14 



